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Abstract 



We introduce a new simple hierarchically constrained model of slow relax- 
ation. The configmational energy has a simple form as there is no coupling 
among the spins defining the system; the associated stationary distribution 
is an equilibrium, Gibbsian one. However, due to the presence of hierarchical 
constraints in the dynamics the system is found to relax to its equilibrium 
distribution in an extremely slow fashion when suddenly cooled from an ini- 
tial temperature, Tq, to a final one Tf. The relaxation curve in that case can 
be fit by an stretched exponential curve. On the other hand the relaxation 
function is found to be exponential when Tf > Tq, with characteristic times 
depending on both Tf and Tq, with characteristic times obeying an Arrhenius 
law. Numerical results as well as some analytical studies are presented. In 
particular we introduce a simple equation that captures the essence of the 
slow relaxation. 
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I. INTRODUCTION 



It has been recognized that glasses show many interesting universal properties which 
have not been satisfactorily explained yet [0,0. It is well known that glassy states are not 
ergodic ; for example, the diffusive motion of molecules in a liquid near the glass transition 
temperature becomes much slower than the experimental time scale, and the state of the 
system is not able to explore the whole phase space within the observation time scale. The 
relaxation function, q(t), that describes the decay of the system energy (or, analogously, 
of any other relevant magnitude) towards its stationary state, is observed to evolve in an 
extremely slow fashion in the vicinity of the glass transition, and is empirically found to be 
well fit by the Kohlrausch- Williams- Watts or stretched-exponential law, i.e. 

q{t) = Aexp[-b*{t/Tf] (1) 

with exponent < /? < 1 [Q. It is often observed that the exponent j3 decreases with 
decreasing temperature [^]. 

Nearby the glassy transition standard methods of statistical mechanics based on equi- 
librium distributions are no longer suitable to describe this kind of systems. A dynamical 
approach is required for an adequate understanding of slow relaxation processes. 

To the date, several dynamical models have been proposed to get some insight into the 
nature of the slow relaxation processes in glassy dynamics. All of them have in common the 
presence of some kind of dynamical frustration, and can be classified in two groups: those 
involving frustration due to the presence of energy and/or entropy barriers (for example spin- 
glass type of models and others in which even with a simple free energy landscape 
the frustration is directly introduced in the dynamics [|10|. The model we present in this 
paper belongs to the second type, i.e. the energy function is very simple, but the dynamics 
is strongly constrained; only a reduced number of degrees of freedom can evolve freely at a 
given time step, while the rest remain frozen. 

The new ingredient we include with respect to previously studied models is the presence 
of a hierarchy of constraints ||TI|]. With that hierarchy we pretend to mimic the fact that 
in slow relaxation processes, usually there are some degrees of freedom that evolve faster 
than others and that influence the dynamics of the slower modes in such a way that a whole 
hierarchy of constrains is generated. 

The idea that a theoretical model of slow relaxation should be a dynamical and hier- 
archically constrained one is not new. It was first introduced by Palmer et. al. |T2|. They 
defined a family of models consisting of discrete levels of degrees of freedom, = 0, 1, 2, . . ., 
each level containing Nn spins. A spin at level n + 1 can change freely its state only when 
A^n < Nn spins in level n happen to be at a given specific configuration among all the 2^^" 
possible ones. As a consequence of such type of dynamical constraints, they argued that a 



stretched-exponential behavior, Eq.(|T]), can be reproduced under broad conditions |T^. One 
important point, however, is that in the previous paper analysis it is implicitly assumed that 
all the spins above a given one in the chain influence the dynamics of it, and therefore the 
constraints are long-ranged. 

In this paper we consider a very simple microscopic hierarchically constrained model for 
slow relaxation, which is a particular physical realization of the general scenario introduced 
with two important differences: 
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1- Our model does not consider long range dynamical constraints, i.e, a degree of freedom 
is constrained directly only by the nearest degrees of freedom in the hierarchy. As we will 
show this is enough to generate a slow relaxation. 

2- The main advantage of this new model is that the dynamics and the constraints are 
specifically defined, rendering the model suitable to be analyzed using computer simulations 
and detailed mathematical analysis. 

The paper is structured as follows: in section 2 we introduce the model, in section 3 we 
present our main numerical results as well as a simple theoretical approximation that repro- 
duces stretched-exponential type of relaxation; we also discuss the nature of the boundary 
conditions. Finally we present the conclusions. A mathematical formulation of the model 
and a mean-field like solution are presented in an appendix. 



II. THE MODEL 

The model consists of a set of degrees of freedom (or spins), evolving according to a 
constrained dynamics. The degrees of freedom could, in principle, be arranged in different 
ways, but in what follows we consider only a one-dimensional chain with a degree at every 
site. Each degree of freedom is exposed to the action of an external potential which functional 
form is schematically shown in Fig.|I|. It has two local minima at positions s = and s = 1 
with values and e, respectively, and a maximum of height cj) > e separating the two minima. 
We assume that each degree of freedom can be located only at the minima of the potential, 
and its state can therefore be characterized by a spin-like variable, s = 0, 1. 

The spins can change stochastically their state with the following (Arrhenius) transition 
rates, V 

To^i = T-o"^exp(-;^) , Ti^o = T-o"^exp(-^—- ^) , (2) 

where Tq is a constant that defines a microscopic time scale, and Tj denotes the temper- 
ature. Assuming that there is no static interaction among spins, we can easily calculate 
the equilibrium probability Peq to find a spin at state s = 1 by using the detailed balance 
condition: 

^ To^i + ri^o ^ exp(e/T^) + 1 

which gives the Gibbs equilibrium distribution associated to the potential in Fig. 1 (observe 
that it does not depend on 0). In the absence of interaction among spins the system exhibits 
exponential (Debye) relaxation towards the equilibrium distribution with characteristic re- 
laxation time given by: r = ppeg^o exp((/)/T/). 

We now introduce a dynamical constraint in the model: in Fig.|] (a) we show schemati- 



cally the way in which a spin, 5*1, constrains the motion of its right neighbor 5*2 [0. When 
the spin 5*1 is at the state Si = 0, the dynamics of spin 5*2 is obstructed, i.e., (p ^ oo, and 
jumps of 5*2 between the two states are completely prohibited. In the case in which both 
spins Si and 5*2 are initially at state s = 0, for 5*2 to jump to the state S2 = 1, 5"! has 
to evolve first to Si = 1 to clear the path of 5*2. Generalizing this rule the motion of the 
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{i + l)-th spin in the chain is blocked by the i-th spin when this one is at state Si = 0. 
Note that while the motion of Sj is affected by Si_i it does not depend on Sj+i, originat- 
ing in this way a hierarchy of directed constraints. It is important to observe that the 
dynamics is asymmetric, and that even though there is no static interaction among spins 
they are dynamically constrained; therefore the equilibrium distribution is the product of 
the single-site equilibrium distributions eq.(§). This is rigorously proved in the appendix 
where a probabilistic formulation of the model is presented. 

Due to this series of constraints we expect the system to show slow relaxation. One 
way to see that is in terms of the topology of the phase space: as discussed before, in the 
two-spin case the system cannot go directly from {si, S2} = {0, 0} to {si, S2} = {0, 1}, but 
has to take a roundabout path {0,0} {1,0} {1, 1} {0, 1}. In the case of L-spin 
system the phase space forms a L-dimensional hypercube with sides of size 1. A vertex of 
the cube represents a system state. Though the maximum geometrical distance between 
any two states is L2 measured by the edge length of the hypercube, a big fraction of the 
possible paths are obstructed and unavailable. Therefore the system has to pass through 
more complex roundabout paths as the system size L increases, and we expect those complex 
trajectories to give raise to slow relaxation. 



III. RESULTS 

We have performed numerical simulations of the model for different system sizes, L, up 



to L = 1000, with open boundary conditions [|Tj], i.e. the first spin in the chain (hierarchy) 



evolves in an unconstrained way (this mimics the fact that typically in glasses there are 
fast, unconstrained, degrees of freedom). Most of the plots we present are obtained for 
L = 250, but the results have been observed to be very robust when increasing system size. 
Simulations have been carried out for both simultaneous and sequential type of update, and 
essentially no physical difference is observed among them. Some parameter values are kept 
fixed in all the simulations: tq = 10.0, (f) = 1.0, and e = 0.5. We have verified that the 
qualitative general features exhibited by the model do not depend on the choice of these 
values. As initial condition the spins are placed at positions s = or 1 with probabilities 
corresponding to an equilibrium distribution, Peg(To) eq.(|), for a given initial temperature, 
Tq. Therefore the dynamics drives the system from an equilibrium distribution at Tg to a 
different equilibrium distribution at Tf. 

We introduce the mean energy per spin defined as 

U{t)^L-hj:s.{t) (4) 
1=1 

(only spins in the s = 1 state give a non- vanishing contribution to the energy). Given 
that there is no configurational interaction among spins, after sufficiently long times, U{t) 
approaches its equilibrium value Ueq = Peq^- We define a relaxation function q{t) as 

q{t)^\U{t)-U,^\ . (5) 

which, after a sufficiently long time, t, q{t) approaches zero. 
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A. General features of the relaxation curves 

Two typical relaxation curves are shown in figure ^ they correspond to a case in which 
the system is cooled down, Tj < Tq (upper one), and in the other the system is heated up, 
i.e., Tf > Tq. Their behaviors are essentially different. 

The uppermost corresponds to I/Tq = 2.19, 1/T/ = 3 and L = 250. Note that the 
relaxation is very slow; as the time is measured in units of Tq = 10, the maximum time 
corresponds to 15000 time Monte Carlo steps (the curve is the average of 10^ independent 
runs). Observe that asymptotically, i.e., after a transient of about tc ~ SOOtq time steps the 
curve is locally well fit by an exponential with a very large characteristic time r ~ 2100ro, 
that implies an extremely slow relaxation. However, it is not clear from numerics whether 
for larger times this exponential behavior will persist, or the curve will decay even in a 
slower fashion. On the other hand, the whole curve, including the initial faster decay can be 
perfectly fit by an stretched exponential with /5 = 0.38 (in fact the fit in undistinguishable 
from the numerical data in Fig. 

Increasing further the time in the computer simulation to decide whether the real asymp- 
totic behavior is an exponential or an stretched exponential is beyond of our available com- 
putational power. In any case, for situations in which the system is cooled down, we always 
get curves that bend progressively in a semilogarithmic plot; and even though the final part 
can always be fit with a straight line (exponential behavior) in no case it is evident whether 
that exponential fit gives the right asymptotic behavior. Even for extremely long times we 
have this type of ambiguity. However, as for any "reasonable" time a stretched exponential 
can always be fit, we admit the relaxation to be non-exponential is this case. At any rate, 
the relaxation is extremely slow in the cooling case. 

The lower curve in figure ^ corresponds to the same initial temperature, I/Tq = 2.19, 
and system size, than the upper one, but a larger final temperature: 1/Tj = 1. In this 
case the transient is much smaller, tc ~ 60, and after it, a pretty clear exponential behavior 
settles in. 

The situation described for the two previous examples is general for all the relaxation 
curves: those in which Tj < Tq are well fit by a stretched exponential, while for the opposite 
situation, Tf > Tq, the decay is exponential after a transient. 

We now analyze the transition between the two previous regimes in a more quantitative 
way. In figure ^ we show the transient time (i.e. the time after which an exponential fit is 
adequate) as a function of the final temperature for a fixed initial temperature 1/Tq = 2.19. 

Note that for very large Tf (small 1/Tj) the system can relax very fast, there is no 
effective frustration and tc is small. In fact, for 1/Tj the system has a huge degree of 
thermal activation, and decays exponentially fast to the equilibrium state with no transient. 
At Tf = Tq the system is already at equilibrium therefore tc = 0. Between these two limiting 
case tc is larger than and behaves in the form shown in Fig. |^; observe that in this interval 
tc is always relatively small, therefore the asymptotic exponential settles in a short time 
(smaller than the observation time). On the other hand, when Tf < Tq (1/7/ > 1/Tq), 
the behavior is quite different and tc grows fastly without bound. We point out that the 
previous set of values of tc for Tf < Tq are obtained fixing some maximum time (t=1500). If 
that time is diminished the apparent tc decreases, while it seems to increase for longer times. 
In this way the values plotted in Fig. ^ are lower bounds for the transient times. In fact, as 
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discussed above, if tc converged to a fixed value for t — oo, that would imply an exponential 
asymptotic decay; while a continuously growing would imply a stretched exponential type 
of asymptotic decay. Deciding which of those two possibilities is the right one from numerics 
is a very difficult task given the extreme slowness of the relaxation process. 

Let us now ignore the transients and present a description of the system behavior in 
terms of its asymptotic decay for the Tf > Tq case. 

For a given pair of fixed initial and final temperatures, the relaxation function is well fit 
(after the transient) by an exponential (as discussed previously) with a characteristic time 
r that depends on both the initial and the final temperature. 

In Fig. ^ and ^ we present the results of our simulations for a chain of length L = 250 
(the results do not change significatively with increasing system size). 

In Fig. ^ the dependence of the characteristic time on the final temperature is shown. 
The curve is well fit by an Arrhenius law: r = Ci exp(— C2/TJ) where Ci and C2 are 
constants. 

In Fig H we keep fixed the final temperature, l/Tj = 2 to study the dependence of 
the characteristic time on the initial state, we also observe that an Arrhenius law: r = 
Di exp(— D2/T0), with Di and D2 constants holds. Therefore the characteristic time depends 
on both the final temperature and the initial state. 

We want to point out that in the two previous graphs, if we crossed to the other regime 
(i.e. Tf < Tq), plotting the characteristic time associated to a long time simulation (for 
example, the one used in figure |^) we would obtain a discontinuous jump: the characteristic 
time for a slight system cooling {Tf = Tq — a, with a small and positive) is much larger (if 
any) than the one for a slight heating up (Tf = Tq + a). 

Physically the main difference among the cooling and the heating processes is the follow- 
ing: as can be easily derived from the equilibrium distribution, the average length, < I >o, 
of chains of spins in the blocking position, s = is: < / >o= 1 + exp(e/T). Observe that it 
goes to 2 at infinite temperature and diverges for vanishing temperature. Analogously the 
mean length of s = 1 chains is < / >i= l + exp(— e/T). In order to cool the system down the 
typical length of a chain of blocking states has to be enlarged. But, of course, the dynamics 
is restricted to sites preceded by islands of s = 1 states. As the number and typical length 
of these islands with s = 1 is decreased in the cooling process the dynamics becomes slower 
an slower. On the other hand, when the system is heated up, there is no such an effect: the 
dynamics is accelerated as more and more active sites appear in the system configuration. 
In order to render more explicit the previous argument let us perform a mean field type of 
calculation: as we have seen the total energy of the system is proportional (through e) to the 
number of spins Ni{t) in the state 1. For A'^i(t) we can write a mean-field master equation: 

m{t + 1) = N,{t) + {L- m{t))r],{t) - m{t)r]2{t) (6) 

where L is the system size, L — Ni{t) is the number of spins in the state at time t, r]i{t) 
is the probability that at time t a spin in the state fiips to 1 and r]2{t) is the probability 
of fiipping from 1 to 0. If we now divide eq.|^ for L, we obtain: 

a{t + 1) = a{t) + (1 - a(t))r/i(t) - a{t)r]2{t) (7) 

where a{t) is the averaged density of spins at the s = 1 state at time t. In this way, eq. 
(0) is nothing but an effective balance equation, which establishes that the probability of 
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being in s = 1 at time t + 1 is given by the probability of being at s = 1 at time t, plus the 
probability of being at s = times the probability of jumping from s = to s = 1, minus 
the probability to escape from s = 1. 

We could consider rji and 772 as given by some fixed values, but in order to keep more 
track of the original constrained nature of the dynamics dynamics we choose them in a 
slightly different way: as a{t) can be interpreted as the probability that a fixed spin is in 
the state s = 1 at time t, we take rji = with probability 1 — (T{t), that is the transition 
probability is if the preceding spin is in a blocking state, and rji = exp(— 1/Tj) (we have 
taken = tq = 1) with complementary probability ait), according to eq. (Q). Analogously, 
?72 = with probability 1 — cr(t), ?72 = exp((e — with probability cr{t). 

In this way, as the system relaxes to its stationary state, the transition probabilities 
change and modify the rate at which the dynamical variable a{t) evolves. 

The previous simple (but not trivial) equation can not be solved analytically, but it is 
simple to iterate numerically. We have considered a{t = 0) given by its equilibrium value 
at some temperature Tq, and iterate the equation for different values of the temperature Tj. 
The result of such an analysis are shown in figure |^. For any set of parameter values the 
resulting relaxation curve can be fit by an stretched exponential with exponent (3 < 1. We 
will further study the curious analitycal properties of this simple approach in a future work. 

Therefore, even in this simple approach, in which the spatial structure is almost com- 
pletely disregarded, we reproduce a stretched-exponential tj^e of decay, due to the fact that 
the jumping probability decays in time as the number of blocking states grows. 



A possible way to synthesize a stretched-exponential is by the convolution of a number 
of exponential curves with different characteristic times, Tn, i.e. 



where Wn are some weight factor. 

In order to shed some light on the microscopic origin of that slow relaxation process, 
and clarify whether the stretched-exponential behavior describing the transient behavior 
of the relaxation function comes from a convolution like that of Eq.(|D (with the n in the 
sum being the position in the chain), we have studied the time relaxation of the energy 
of the spin at every site as a function of its position in the chain. In Fig.^ we show the 
time evolution of U{i,t) — Ueg, with U{i,t) = e{si{t)), where (. . .) stands for averages over 
different runs. It is observed that the first spin in the chain relaxes faster than the second 
one, the second one faster than the third, and so on. I fact, the energies of the first spins 
U{i, t) relax exponentially fast to its equilibrium value, Ueq = Peq^, with a time constant, 
Ti, that increases with the position i in the chain. After a certain number of spins that 
depends on the parameter values, the relaxation of every single spin is indistinguishable 
from the relaxation observed performing a simulation with periodic boundary conditions 
and the same parameter values: this is what we call the bulk behavior. Therefore, the only 
difference between periodic and open boundary conditions is a small effect that does not 
affect the bulk properties. As apparent stretched-exponential behavior are observed for the 



B. Open versus periodic boundary conditions 




(8) 
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relaxation of the bulk spins when Tf < Tq, wc conchide that it is not due to a convohition of 
exponential functions with different characteristic times associated to the different positions 
in the chain. 

An alternative possible way to understand the apparent stretched exponential behavior as 
associated to a convolution of exponentials with different characteristic times is by assuming 
that the islands of spins in position s = 1 in the system relax in a different way depending 
on how many blocking spins are placed in the immediately superior positions in the chain. 
In this way it is clear that spins blocked by only one spin will relax much faster than spins 
preceded by long chains of blocking spins. This mechanism would give rise to different 
relaxation times for spins in different relative positions in the chain, and therefore to a 
global stretched-exponential behavior. 

IV. CONCLUSIONS 

We have presented a simple model of slow relaxation. Its stationary equilibrium distri- 
bution is a simple Gibbs distribution, but its dynamics is strongly constrained. That gives 
rise to slow relaxation processes in the case in which the system is cooled from an initial 
temperature, Tq to a smaller Tf. The asymptotic behavior of the relaxation function, after a 
transient, is found to be exponential when heating the system up, with a relaxation time de- 
pends on both the initial and the final temperatures in a non-trivial way. When cooling the 
system down it is not clear whether the relaxation function reaches an exponential behavior 
asymptotically or if it is a stretched-exponential even asymptotically. In any case, for small 
enough final temperature, and for large initial temperatures, the system is found to show a 
extremely slow relaxation originated by the constrained dynamical rules. In this way we put 
forward, using a very simple model, how constrained dynamics can slow down relaxation 
processes in quite a dramatic way. We have also introduced a very simple one- variable equa- 
tion that captures the essence of the slow relaxation and reproduces stretched-exponential 
type of decay. 

V. ACKNOLEDGEMENTS 

We thanks C. A. Angell, U. Marini BettoUo, Sebastiano Carpi, and G. Parisi for useful 
discussions. We are grateful to the anonymous referee whose comments and criticisms helped 
us to improve notoriously this paper. This work was partially supported by the Japan Society 
for the Promotion of Science through a grant to H.I. and by the European Union through a 
grant to M.A.M. 

VI. APPENDIX: ANALYTICAL APPROACH 

Now we present an analytical attempt to understand the previously described properties. 
Our model can be formally represented in terms of a Markov chain, defined by the following 
equation: 

P{{sh = E WiW} - {s})P{{s'}; t - 1) (9) 

{s'} 



8 



where {s} = {si = Si{t)}, and {s'} = {s'i = Si{t — 1)} for every i = 1,2, ...,L, are system 
configurations at time t and t — 1 respectively. The equation ^ states how the probabihty, 
P{{s}; t) of finding the system at a given time t in a configuration {s} evolves in time. The 
transition probabilities, {s})P({s'};t — 1) are given by: 

L 

W{{s'} {s}) = n u^is', ^ s,; (10) 

i=l 

with 

roexp(/?0) 

and the boundary condition So(t) = 1. This is nothing but the mathematical expression of 
the transition probabilities described in the previous section. 

By direct substitution it is easy to verify that the equilibrium distribution 

is the stationary solution of Eq.(|D. 

From the general equation Eq.(|^), we can derive a hierarchy of equations for the m- 
body probability distributions (similar to the BBGKY hierarchy in statistical mechanics). 
In particular, for m = 1 

1 1 

p{si]t) = Yj 'Y uj{s'i ^ Si] s'i-i)p{s'i-i,s'i]t - 1) (13) 

s'i_i=Os'i=0 

where 

p{sf,t) = j: ... j: E ... EP(M;i) (14) 

si=0 Si_i=Osi+i=0 sl=0 

are the one-body probability functions at time t, 

p{s',.us\;t-i)= j: ... E E ... EP(M;^-i) (15) 

s'i=0 s',_2=0s'i+i=0 s'l=0 

are the two-body probability functions at time t — 1. Analyzing this set of equations is a 
difficult task as can be seen from the fact that the one-body probability functions depend 
on the two-body probabilities distributions p(si_i, sf, t — 1); the equations for these depend 
on three-body probabilities and so on. 

The first analytical approach ( the only one we complete in this paper) consists of approx- 
imating the two-body probabilities by the product of two one-body probability functions, 
namely: 

p{s'i_us'f,t-l)^p{s'i_^;t-l)p{s'i;t-l) . (16) 
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This is a mean field like approximation given that high order correlation are neglected. 
Introducing Eg. ([TBI) in Eg . (|13|) , and after some simple algebra, we get 



P^{t) -Pi{t - 1) = m-^[P^-l{^ - - 1)) - exp(/3e)pi_i(t - l)p,{t - 1)] (17) 

where pi{t) = p{si = 1; t). In the case of an open chain we put poit) = 1, which corresponds 
to the fact that the first spin is unconstrained at any time. 

In the continuous time limit, i.e. when tq ^ 1 (which is the case in the numerical 
simulations), the previous eguation can be written as: 

^ = ^^^^p(^^) [P.-i(t)(l -P.(t)) -exp(/je)p.^,(t)p.(t)] (18) 

which general solution is 

Pi{t) = Peq + (Pi(0) - Peg) exp(/j (t) /Ti ) , (19) 

with Ti = TQPeqexp((f)/Tf ) (which is the free relaxation time) and Jj = Jq Pi-i{t') dt' is an 
effective time for i-th spin. Substituting Ji = t in Eg . (|T9|) , we can calculate J2. It is obtained 
to be given asymptotically by I2 — Peqt for t ^ T/p^q, and therefore Eg. (p!9|) , for i = 2 can 
be rewritten as 



Piii) = Peq + (Pi(0) - Peg) exp{-pegt / Ti) . (20) 

By iterating the same procedure we get Jj ~ pegt and Eg. (pO|) as the solution for pi for any 
i > 1. The relaxation time for i-th spin Xj by Eg.(^) is 

T"i = Ti/peq ~ exp(const. X ^) , (21) 

which has a temperature dependence as that of strong liguids. However, this mean field 
approximation does not reproduce correctly the characteristic decay times. We expect that 
more accurate approximations to Eg.(^ reproduce more accurately the behavior of the 
relaxation function. Results obtained by truncating the hierarchy at higher levels will be 
presented elsewhere [ pTS] . 
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FIGURES 
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FIG. 1. A schematic figure of the potential form in which atoms are placed in our model. Two 
local minima s = and 1 are separated by a maximum of height (f). 

FIG. 2. When the spin 1 is placed at position si = the transition of the spin 2 is prohibited. 

FIG. 3. The phase space for the two spin system. The direct path from {si,S2} = {0,0} to 
{0, 1} is obstructed and the transition between these two states is prohibited. 
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FIG. 4. Typical relaxation curves for Tf < Tq (uppermost, l/Tj = 3 and I/Tq = 2.19), and 
Tf > To (lowermost 1/7/ = 1 and 1/Tq = 2.19). The second one is exponential while the first one 
can be fit by an stretched exponential with /? = 0.38; an exponential fit (dashed line) is valid for 
times larger than ^ 800. L = 250. 




FIG. 5. Semilog plot of the inverse of the characteristic time as a function of l/Tj for a fixed 
l/To = 2.19. L = 250. 
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FIG. 6. Semilog plot of the inverse of the characteristic time as a function of I/Tq for a fixed 
l/T/ = 2. L = 250. 
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FIG. 7. Transient time (as defined in the text) as a function of 1/Tf, with I/Tq = 2.19. 
L = 250. 
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FIG. 8. Stretched exponential decay associated to eq.(7), when the system relaxes with 
l/Tf = 3 from an initial (t(0) = 1/(1 + exp{l/2)). The fit is given by 6 = 3.7535 and (3 = 0.201 
(see eq.(l)). 
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FIG. 9. Semilog plot of the decay of f/(i, t) — Ueq for the first seven spins in an open chain as 
a function of time, for a system with L = 250. The lowermost curve corresponds to the first spin 
in the chain, the next one to the second one and so on. 
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